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ABSTRACT 

Primordial star formation appears to result in stars at least an order of mag- 
nitude more massive than modern star formation. It has been proposed that the 
transition from primordial to modern initial mass functions occurs due to the 
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onset of effective metal line cooling at a metallicity Z/Z = 10 -3,5 . However, 
these simulations neglected molecular hydrogen cooling. We perform simulations 
using the same initial conditions, but including molecular cooling, using a com- 
plex network that follows molecular hydrogen formation and also directly follows 
carbon monoxide and water. We find that molecular hydrogen cooling allows 
roughly equivalent fragmentation to proceed even at zero metallicity for these 
initial conditions. The apparent transition just represents the point where metal 
line cooling becomes more important than molecular cooling. In all cases, the 
fragments are massive enough to be consistent with models of primordial stellar 
masses, suggesting that the transition to the modern initial mass function may 
be determined by other physics such as dust formation. We conclude that such 
additional cooling mechanisms, combined with the exact initial conditions pro- 
duced by cosmological collapse are likely more important than metal line cooling 
in determining the initial mass function, and thus that there is unlikely to be a 
sharp transition in the initial mass function at Z/Z = 10~ 35 . 

Subject headings: stars: formation - stars: mass function - early universe - 
hydrodynamics - equation of state - methods: numerical 



Introduction 



The formation of stars is a key process in the early universe with far-reaching conse- 
quen c es for cosmic r eionization and galaxy formation ( Loeb &: BarkanalboOl : Bromm Sz Larson 
2004 : iGloverl 120051 ). The physical processes that govern stellar birth in a metal-free or 



metal-poor environment, however, are still very poorly understood. Numerical simulations 
of the thermal and dynamical evolution of gas in primordial protogalactic halos indicate 
that the metal-free first stars, the so called Pop ulation HI, are expected to be very massive 



with masses anywhere in the range 20-2000 (I Abel. Bryan, fc Normanll2002l ; iBromm et al. 



20021 : lYoshida et a/.lbOOGuO'Shea fc Normanj|2007l ). In contrast, the stellar mass spectr um in 



Kroupa 


2002: 


Chabrier 


2003) 



there must have been a transition from primordial, high-mass star formation to the "nor- 
mal" mode of star formation that dominates today. The discovery of extremely metal-poor 



subgiant stars in the Gal actic halo with masses below one solar mass (IChristlieb et al.l 12002 



Beers fc Christliebl 120051 ) indicates that this transition occurs at abundances considerably 
smaller than the solar value. At the extreme end, these stars have iron abundances less 
than 10" 5 t imes the solar value, and carbon or oxygen abundances that are still ^10 
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times the solar value. These stars are thus strongly iron deficient, which could be due to 
unusual abundance patterns produced by substantial mixing a nd fallback of ejecta dur ing 
supernovae explosions of stars with masses between 2 0-130 M m dUmeda fc Nomotoll 20031) or 



due to mass transfer from a close binary companion (IRyan et al 



2005 



Komiva et all 120071 ). 



There are hints for an increasing binary fraction with decreasing metallicity for these stars 
jLucatello et alil2005l ). 



In a seminal paper, iBromm et al.l (120011 ) proposed that the transition from high to 
low mass star formation should occur at a critical metallicity Z cr ;t ~ 10 ■ Zq as a result of 
atomic fine-structure line cooling from alpha elements, such as carbon and oxygen, becoming 
important. Their study is based on smoothed particle hydrodynamics simulations of the 
collapse of cold gas in a top-hat potential at densities below n = 10 6 cm~ 3 . It assumes that 
cooling from molecular hydrogen is negligible in the density and metallicity range considered. 
The value of Z cr - lt can also be estimated analytically by calculating the metallicity required 
to produce a cooling rate equal t o the rate of adiabatic c ompression heating for Riven hal o 



prop ertie s (IBromm fc 



ft 

oebl 120031 : ISantoro fc Shulll 120061 : iFrebel. Johnson fc Bromml I2007D . 



However, IBromm et al.l (l200ll ) noted that neglecting H2 cooli ng could signifi c antly influence 
the resulting fragmentation pattern. Along the same lines, lOmukai et al.l (120051 ) argued, 
based on the results of detailed one-zone calculations, that molecular cooling would indeed 
dominate the cooling over many orders of magnitude in density. They argued that dust- 
induced cooling at densities above n « 10 13 cm -3 may be a more viable mechanism to explain 
fragmentati o n in y ery low metallicity gas, a conjecture that has recently been confirmed by 
Clark et all (I2008I ). 



The effect s of m olecular cooling at densities up to n rs 500 cm -3 have been discussed by 
Jappsen et al.l (120071 . hereafter Paper II) in three-dimensional collapse simulations of warm 
ionized gas in protogalactic halos for a wide range of environmental conditions. This study 
used a time- dependent chemica l netw ork running alongside the hydrodynamic evolution as 
described in iGlover fc Jappsen! (120071 . hereafter Paper I). The physical motivation was to 
investigate whether small protogalaxies that formed within the relic Hn regions left by neigh- 
boring Population III stars could form subsequent generations of stars themselves, or whether 
the elevated temperatures and fractional ionizations found in these regions suppressed star 
formation until larger protogalaxies formed. We found in Paper II that molecular hydrogen 
dominates the cooling of the gas for abundances up to at least 10~ 2 Z . In addition, there 
was no evidence for fragmentation at densities below 500 cm -3 . 



To understand whether this apparent discrepancy with the IBromm et al.l (1200 ll ) re- 
sults comes from including molecular hydrogen cooling or from adopting vastly different 
initial conditions — here ionized halo gas in a dark matter halo with a NFW density profile 
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( INavarro et al.lll997l ) , there cold atomic gas in a top-ha t halo — we apply ou r time-dependent 



chemical network with molecular hydrogen cooling to iBromm et al.l (l200ll )-type initial con- 
ditions. Our study is structured as follows. In the next section, §2], we give the details on 
the numerical set-up and the adopted initial conditions which we summarize in Table [TJ In 
£J3] we report our results, and we summarize in §H 



2. Numerical Set-Up 



To adequately describe star formation in the early universe, it is necessary to follow 
the cooling and possible fragmentation of the gas in the central regions of dark-matter halos 
over many orders of magnitude in density. Due to the nonlinear nature of the flow, it is not 
known in advance where and when gravitational collapse will occur, and consequently stars 
will form. To compute the time evolution of the system we therefore resort to smoothed 
particle hydro dynamics ( S PH), which is a Lagran gian scheme to solve the equations of hy- 
drodynamics (Benz 1990 ; Monaghan 1992^ 2005 ) coupl ed with the a simplified version of 
the time-dependent network discussed in iGloverl (120091 . hereafter, Paper III). The fluid is 
represented by an ensemble of particles, and flow quantities are obtained by averaging over 
an appropriate subset of SPH particles. The method is able to resolve high density contrasts 
as particles are free to move, and so the particle concentration increases naturally in high- 
densi ty regions. We use a modified version of the parallel SPH code GADGET (jSpringel 



20051 ). Once the central pa rts of a collapsing fragment exceed a density n = 10 5 cm 3 we 



introduce a "sink particle" (IBate et al.lll995c IJappsen et al.ll2005l ). which is able to accrete 
gas from its surroundings while keeping track of the mass and linear and angular momentum 
of the infalling material. Replacing collapsing high-density cores by accreting sink particles 
allows us to follow the dynamic evolution of the system over many local free-fall timescales. 



2.1. Chemistry and Cooling 

The chemical model used in our simulations is based on the detailed model of low- 
metallicity gas chemistry outlined in Paper III. This model includes cooling both from atomic 
fine structures lines and also from the molecules H 2 , HD, CO, OH and H 2 0. It was designed 
to accurately follow the formation and destruction of these molecular species over a wide 
range of temperatures and densities. 

In Paper I, we showed that the neglect of molecular species including metals can be 
justified in simulations of gravitationally collapsing protogalactic gas provided that n < 
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500/t cri ar cm -3 , where where t c har is the characteristic physical timescale of interest, in Myr. 
At higher gas densities, molecule formation from metals may become more significant, al- 
though its ultimate importance depends on a number of factors, such as the temperature 
and metallicity of the gas, and the presence (or absence) of an ultraviolet radiation field. 

Although our chemical model is capable of accounting for the effects of dust, cosmic 
rays, and photodissociation by an external radiation field, we omit these complications in 
the present study. 

A full list of the chemical reactions included in the chemical network used in this paper 
is given in Tabled along with the corresponding rate coefficients and references to the source 
of the data. Given these rate coefficients, calculation of the appropriate chemical rates is for 
the most part straightforward. One minor exception is the case of H 2 collisional dissociation, 
which has a rate coefficient that depends on both temperature and density. For the collisional 
dissociation of H 2 by H or H 2 (reactions 9 & 10), we calculate the rate coefficients using a 
function of the form 



log h 



n/n c 



log fc if h + 



log h t ] 



V. 



l + n/n CI J ' \l + n/n. 

where k- h \ and k^ are the low density and high density limits for the collisional dissociation 
rate due to collisions with species i (which can be found in Table [2]), and where the critical 
density, n cr , is given by 

(2) 



Xh 2 



™cr,H n cr,H 2 

where n CTj H and ri C r,H 2 are the critical densities in pure atomic gas with an infinitesimally 
dilute quantity of H 2 and in pure molecular gas respectively. Temperature dependent fits 
for n cr H and n cr) H 2 can be found in Table [2J For the collisio nal dissociat i on of H 2 by He 
(reaction 75), we use a slightly different expression, taken from lDove et al.l (119871 ): 

(log /c He ,h - log /c He ,i) 



log k Re = log fc He ,h 



\1.09 



(3) 



1.0 + (n/n CTjHe 

Values for fcHe,i, fcHe,h and n cr He are given in Table El 

A referenced list of the thermal processes included in our models is given in Table [3j 
To compute the CO and H 2 rot ational and vibrational co oling; rates, we use the multi- 
param eter fitting functions given in lNeufeld fc Kaufman! (119931 ) and lNeufeld. Lepp &: Melnick 
(119951 ) . One of the parameters required is the effective column density N(M) of each coolant 
M (where M = CO or H 2 Q). We calculate this in our simulations using the expression 



(INeufeld & Kaufmanlll993[ ) 



N{M) 



n(M) 
IV- v| 



(4) 
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whe re n(M) is the local num b er de nsity of species M. The iNeufeld fc Kaufman! (119931 ) 
and INeufeld. Lepp fc Melnick! (119951 ) treatments assume that only collisions with H 2 are 
important in determining the CO or H 2 rotational cooling rates, as is appropriate in a 
fully molecular gas. However, in metal-poor gas, the H 2 abundance will often be much 
smaller than the abundances of atomi c hydrogen or atomic helium. We therefor e use an 
effective number density n e s, given by (IPavlovski et al.l 120021 ; ISmith fc Rosen! 120031 ) 



n c g = nun + V^nu + 0.5n 



He 



(5) 



i n plac e of the H 2 number density in the INeufeld fc Kaufman! (119931 ) and lNeufeld. Lepp fc Melnick 
( 119951 ) fitting formulae. 



Radiative cooling alone is unable to decrease the temperature of the gas below the 
cosmic microwave background (CMB) temperature, on simple thermodynamical grounds. 
We account for this effect in an approximate fashion by using a modified radiative cooling 
rate in our treatment of the gas thermal energy equation. This modified rate, A mod (T), is 
given by 

A mod (T)=A(T)-A(T CM B), (6) 

where A(T) and A(Tcmb) are the unmodified radiative cooling rates at T and Tcmb respec- 
tively. 

For additional details of our chemical networks and cooling function, please consult 
Papers I and H I; further deta i ls of their implementation within Gadget are also given in 
Paper II and in IJappsen et al.l (120081 . hereafter Paper V). 



2.2. Initial Conditions 



We model the initial dark matter distribution, following iBromm et al.l (120011 ). as an 
isolated and roughly spherical overdensity described by the top-hat approximation with 
small-scale fluctuations with an initial variance of a 2 ~ 0.1 added. We consider a halo 
virializing at z w - 1T = 30, which corresponds to a roughly 3cr peak with a total mass of 2 x 
10 6 M Q , corresponding to 2 x 10 5 M Q in baryons. The overdense region has a proper radius 
of 150 pc. For mor e details on the set-up of the dark matter and the gas distribution see 
Bromm et all (boOlh . 



We carry out numerical experiments with five different metallicities ranging from Z = 
to Z = 0.1 Z . Aside from the difference in metallicity, these calculations have identical 
initial conditions. The simulations were initialized at a redshift z = 100, with an initial 
gas temperature T gaSj i = 200 K. To study the effect of resolution we ran simulations with 
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three different resolutions. Our low resolution simulations, have N~dm = 13400 dark matter 
particles and N gas = 65536 gas par ticles. In thes e simu lations, our numerical resolution of 
M res = 150 M was the same as in iBromm et al.l (120011 ). Our higher resolution simulations 
with N BM = 107177 and N gas = 8 x 10 5 have mass resolution M res = 12.5 M©. Finally, 



our highest resolution simulation has ^dm = 4057344 and N gas = 3.2 x 10 6 , so its mass 
resolution is M res = 3.1 M . In all of our simulations, both the dark matter and the gas 
particles were endowed with the sa me Hubble expansion and were set in rigid rotation with 



a spin parameter of 0.05, just as in IBromm et al. 



(1200 if ). 



We also consider a halo viralizing at a lower redshift of z v ir = 25, with a total mass of 
2 x 10 6 M Q , where the overdense region has a proper radius of 200 pc. We summarize the 
properties of all of the runs in Table [TJ where the runs with halos that collapse at a lower 
redshift are denoted with "-C". 



3. Results 

In each of our calculations the halo contracts rapidly and the gas density increases by 
several orders of magnitude. A disk-like structure builds up in the very center of each halo 
with a density of n ~ 10 3 cm -3 and a diameter of 10 to 20 pc, as illustrated in Figured! The 
disk mass at the time illustrated is of order ~ 2-3 x 10 4 M Q , and so it is well-resolved 

even in the low resolution simulations. The disk is supported by a combination of thermal 
pressure and rotation and its growth is fed by material falling in along filaments and sheets. 
This complex disk provides the background for the small-scale initial density fluctuations to 
grow, some of w hich become gravitati onally unstable and collapse to form stars. In contrast 



to the results of IBromm et al.l (120011 ). we find that fragmentation of the disk occurs in all 
of our simulations, and that there is no evidence for a sharp transition in the initial mass 
function at Z = 10~ 3 - 5 Z Q . 

This result can be understood as follows. Regardless of the metal content, H 2 line cooling 
maintains the gas temperature at densities above 10 3 cm -3 close to the CMB temperature 
(2cmb = 105 K at the redshift z = 37.5 considered in Figure [2]); even in the primordial case, 
the gas temperature is within a factor of two of this value. The resulting equation of state 
is therefore approximately isothermal (at high densities) or softer than isothermal (at low 
densities) . The thermal behavior of the gas is depicted in the left column of Figure [2J 

The physical conditions where the equation of state changes from softer than isothermal 
to isothermal or stiffer, i.e. from a regime where the gas cools with increasing density to 
one where the temperature remains constant or rises again, imprint a characteristic mass in 
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the fragmentation process (IJappsen et al.ll2005l ; lLarsonll2005l ). Because this fragmentation 



occurs at approximately the same density and temperature in each simulation, the mass 
distribution of the fragments does not vary greatly from simulation to simulation. In each 
of our low resolution simulations, between 16 and 18 fragments are formed, with masses in 
the range 10 2 M Q < M < 5 x 10 3 M Q . The resulting mass spectra are shown as the solid 
lines in the right-hand column of Figure [2j For comparison, the thermal Jeans mass for gas 
with temperature T gas = T CM b is approximately 2.3 x 10 3 M at n = 10 3 cm -3 , the density 
above which H 2 line cooling maintains the gas as approximately isothermal, and 230 M Q at 
n = 10 5 cm -3 , the density above which we form sinks. To quantify whether the mass spectra 
differ significantly, we can use the two-sided Kolmogorov-Smirnov test to determine for each 
of our low resolution, metal-enriched simulations the probability that the mass spectrum 
produced is drawn from the same underlying distribution as in the run Z0. We show the 
corresponding cumulative distribution functions of the clump masses in the top panel of 
Figure [3J We find probabilities of 0.63, 0.22, 0.99 and 0.11 for runs Z-4, Z-3, Z-2 and Z-l 
respectively, indicating that the differences between the fragment mass spectra are in general 
not statistically significant. We also determine the probabilities for the runs Z-3, Z-2 and 
Z-l that the mass spectrum is drawn from the same underlying distribution as in the Z-4 
simulation. The probabilities are 0.73, 0.63, 0.5, respectively. This shows that there is no 
significant change in the clump mass function at Z = 10~ 3,5 Z . 

In our high resolution simulations, we find a greater degree of fragmentation and a broad 
fragment mass distribution. However, just as in our low resolution simulations, the mass 
distribution is dominated by fragments with masses M > 100 M Q . We find very few objects 
with masses close to our resolution limit of 12.5 M , and there is no indication from our 
current results that we are missing a substantial population of objects with masses below 
this limit. As before, we can compare the fragment mass distributions resulting from our 
three high- resolution simulations with the K-S test. We show our results in the bottom panel 
of Figure [3J We find that the probabilities that the fragment mass distributions in runs Z- 
3-BIG and Z-l-BIG are drawn from the same underlying distribution as in run Z0-BIG are 
0.11 and 0.02 respectively, indicating that there is a statistically significant difference. If 
we compare the runs Z-3-BIG and Z-l-BIG then we determine a probability of 0.93 that 
their drawn from the same underlying distribution. The run with zero metallicity results 
in fragments with masses above 100 M , whereas the runs with metallicity also produce 
fragments with masses between 50 and 100 M . Although we see some slight influence 
from the metal-line cooling, we see no evidence for any low-mass star formation: all of the 
clumps formed are well-resolved, and have masses consisten t with recent determinati ons of 



the plausible mass range for Population III stars (see e.g. lO'Shea fc Norman! 120071 ). We 



see no indication that metal-line cooling at these temperatures and densities is capable of 
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producing low-mass stars. 

We note that the maximum gas temperature during collapse depends on the metallic- 
ity. In the zero-metallicity case, the peak temperature T pea k reached due to compressional 
heating during the initial contraction of the gas is ~ 6 x 10 3 K. As the metallicity increases, 
however, T pea k decreases. In particular, the run with Z = 0.1 Z shows a much smaller peak 
temperature of ~ 10 3 K. This decrease in T peak reflects the growing importance of metal-line 
cooling in cold gas as the metallicity increases. From Figure H] it is apparent that the metals 
begin to affec t T pea \r once the met allicity exceeds Z = 10~ 3 Z . This is a lower threshold 



than found in lJappsen et al.l (120071 ). but the difference is easily understood as a consequence 
of the fact that in these simulations we start with cold gas, and that the H 2 cooling rate 
of the cold, low-density gas in th ese simulations is mu ch lower than H 2 cooling rate in the 



warm, low-density gas studied inlJappsen et al.l (120071) . Our threshold of Z ~ 10 3 Z Q is in 
reasonable agreement with the lBromm et al.l (120011 ) threshold of Z ~ 10~ 3 ' 5 Z ; however, the 



key difference here is that at lower metallicities, the H2 alone is sufficient to cool the gas, and 
the outcome of the fragmentation process has very little dependence on whether the cooling 
comes from the H2 or from the metals. 

We also test the hypothesis that the fragment mass spectra of run Z-l-BIG and of the 
higher resolution run Z-1-BIG2 are drawn from the same distribution and find a probability 
of 0.98. Since the number of fragments are also quite similar we conclude that our runs 
Z-l-BIG, Z-3-BIG and Z0-BIG are converged. 

In all of our simulations the metal-enriched gas is able to cool to the CMB temperature 
of approximately 110 K at this redshift (Figure [2]). This is also true to within a factor of 
2 for the primordial case. H2 can act as a coolant for the gas down to a temperature of 
approximately 200 K, which is close to the CMB temperature at this redshift. Since this 
could be seen as a reason for the lack of influence of the metals that we find, we perform 
two runs Z0-BIG-C and Z-l-BIG- C with halos that collapse at a lower redshift z = 25. In 
Figure H] we show the results of these two runs. At a redshift z = 25 the CMB temperature 
is approximately 71 K and the run Z-1-BIG-C is able to cool down to the CMB temperature 
due to the combination of efficient metal-line cooling and H 2 cooling. In the primordial case 
however the gas is not able to cool at all at this redshift. However, by redshift z = 18, 
a substantial amount of the gas is nevertheless able to cool, collapse, and fragment from 
H2 cooling, even though H 2 cannot cool much below ~ 200 K, and so cannot cool down 
to the CMB temperature at this redshift. We compare the fragment mass distributions 
resulting from the primordial and the metal-enriched case with the K-S test. We find that 
the probability that the fragment mass distributions in runs Z0-BIG-C and Z-1-BIG-C are 
drawn from the same underlying distribution is 0.11. The probability that runs Z-1-BIG-C 
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and run Z-l-BIG are drawn from the same underlying distribution is 0.05 (see Figure [3]). 
The fragments of run Z-1-BIG-C have masses down to 30 M Q whereas the fragments created 
in run Z-l-BIG have masses above 50 M Q . In both runs the gas is able to cool down to 
the CMB temperature due to the combination of metal-line cooling and H 2 cooling. Since 
the CMB temperature is smaller for the run Z-1-BIG-C, the gas can fragment to smaller 
masses. This demonstrates that in runs with higher metallicity, where the gas is able to cool 
down to the CMB temper ature, this tempe rature determines the dynamical behavior (c.f. 
the similar result found by lSmith et al.ll2008l ). Nevertheless, in all of the runs, the fragments 
have masses that are consistent with the expected masses of Population III stars, and there 
is no evidence for the formation of low-mass Population II stars. 



4. Summary and Conclusion 



The current study shows that gas in simulations of the type considered here, with low 
initial temperature, moderate initial rotation, and a top-hat dark-matter overdensity, will 
readily fragment into multiple objects with masses above 30 M , regardless of metallicity, 
provided that enough H2 is present to cool the gas. Rotation leads to the build-up of 
massive disk-like structures in these simulations which allow smaller-scale fluctuations to 
grow and become gravitationally unstable to form protostars. The resulting mass spectrum 
of fragments peaks at a few hundred solar masses, roughly corresponding to the thermal 
Jeans mass in the disk-like structure. 

Our major con clusion is that the metallicity threshold at Z = 10~ 3 ' 5 Z Q reported by 
Bromm et al.l (120011 ) does not represent a critical metallicity above which gas fragments and 
below which it does not. Rather, this threshold simply marks the point at which metal- 
line cooling becomes more important than H2 cooling at the gas densities and temperatures 
relevant for fragmentation in the disk. This metallicity threshold only represents a criti- 
cal metallicity for frag mentation in these sim ulations if H 2 formation is strongly suppressed 
or not considered (e.g. iBromm et al.l 120011). This sce nario may be relevant if the extra- 



galactic UV background is strong (IHaiman et al.ll2000l ). although recent work suggests that 
even in this case, suppression of H 2 formati on is less important than previously supposed 
rtWise fc Abell2007l : lO'Shea fc Normanll2008h . 



We also find that the mass spectrum of fragments formed in these simulations has very 
little dependence on the metallicity of the gas. In our low-resolution simulations, many of 
the fragments are only marginally resolved, but we obtain very similar results from our high- 
resolution simulations, in which the fragments are well-resolved, suggesting that this lack of 
influence is real. 
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We further note that simulations performed by ISmith &: Sigurdssonl (120071 ) that follow 
the evolution of cold, metal-enriched gas from cosmological initial conditions until collapse 
do find tentative evidence for a metallicity threshold at around Z = 10 -3 Z that must be 
surmounted before fragmentation will occur. On the other hand, the simulations presented 
in Paper V that follow the collapse of hot, initially ionized gas into isolated NFW halos 
find no evidence for any fragmentation at or above this metallicity. We conclude that the 
question of whether there is a critical metallicity below which the formation of low-mass stars 
is impossible remains unresolved. Our current results show no evidence for such a critical 
metallicity: fragmentation occurs even in our Z = simulation, due only to efficient H2 
cooling, and the mass spectrum of the fragments formed does not appear to differ significantly 
from those obtained from our higher metallicity simulations. 

However, the degree of fragmentation occurring in these simulations appears to be a con- 
sequence of the initial conditions c hosen: simulations u s ing substantially different initial con- 
ditions find very different results (IJappsen et al.ll2008l ; ISmith fc Sigurdssonl 120071) . Another 
examp le arises if we compare our zero metallicity results with the results of I Abel. Bryan, fc Norman 
(120021 ). These authors use proper cosmological initial conditions and in their simulation only 
one core forms. Thus the fragmentation in our simulation with primordial gas seems to be 
a generic outcome of the adopted initial conditions. It is therefore important to consider 
whether a metallicity threshold appears in simulations with different initial conditions and 
what the most realistic initial conditions are. To make further progress in understanding the 
role (if any) that metal-line cooling plays in promoting fragmentation, we need to develop a 
much better understanding of how metals are initially dispersed into the high-redshift inter- 
stellar and intergalactic medium, in order to be able to select the most appropriate initial 
conditions for our simulations. 



Our results cal l attention to a n alte rnative scenario for the transition to the present- 
day low-mass IMF. lOmukai et al.l (120051 ) stress the importance of dust-induced cooling at 
high densities, which may become im portant at metallicities as low as Z « 10 -6 - 10 -5 Z 
de pending on the adop ted dust model (ISchneider et al.ll2006l ). Recent numerical simulations 
by IClark et al.l (120081 ) support this point of view and predict the existence of low-mass 
stars and even brown dwarfs with metal abundances of ~ 10 -5 the solar value (see also 



Tsuribe fc Omukaill2006h 
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Fig. 1. — From top to bottom: Projected gas number density in runs ZO, Z-3, and Z-l at 
z = 35. Left Face-on view. Right: Edge-on view. The size of the box is 60 pc. 
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Fig. 2. — From top to bottom: Gas properties and mass distribution of clumps at z — 37.5 for 
runs Z0, Z-4, Z-3, Z-2, and Z-l, arranged in order of increasing metallicity. Left: Temperature 
versus number density of the gas. Right: Mass distribution of the clumps represented by 
sink particles. In the panels for runs Z0, Z-3, and Z-l we also show the results for the higher 
resolution runs ZO-BIG, Z-3-BIG, and Z-l-BIG {dotted lines). The dot-dashed line in the 
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Fig. 3. — Top panel: Cumulative distribution function of the clump masses for the low 
resolution runs ZO (solid line), Z-4 (dotted line), Z-3 (dashed line), Z-2 (dot-dashed line), 
and Z-l (dot- dot- dashed line). Bottom panel: Cumulative distribution function of the clump 
masses for the higher resolution runs ZO-BIG (solid line), Z-3-BIG (dotted line), Z-l-BIG 
(dashed line) and Z-1-BIG-C (dot-dashed line). 
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Fig. 4. — Top panel: Gas properties and mass distribution of clumps at z — 18 for the 
run ZO-BIG-C. In this run no clumps have formed yet at z = 25 when the halo starts to 
collapse. Bottom panel: Gas properties and mass distribution of clumps at z = 25 for the 
run Z-1-BIG-C. The dashed line shows the CMB temperature. 
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Table 1. List of the runs discussed in this paper 



Run 


(*©) 




M c 

1 'gas 


-^clump 


ZO 


0.0 


0.05 


65536 


16 


Z-4 


10- 4 


0.05 


65536 


16 


Z-3 


10- 3 


0.05 


65536 


17 


Z-2 


10- 2 


0.05 


65536 


18 


Z-l 


10- 1 


0.05 


65536 


17 


ZO-BIG 


0.0 


0.05 


8 x 10 5 


29 


Z-3-BIG 


10- 3 


0.05 


8 x 10 5 


33 


Z-l-BIG 


10- 1 


0.05 


8 x 10 5 


38 


Z-1-BIG2 


10- 1 


0.05 


3.2 x 10 6 


38 


ZO-BIG-C 


0.0 


0.05 


8 x 10 5 


32 


Z-1-BIG-C 


10- 1 


0.05 


8 x 10 5 


40 



a Metallicity of the gas. 

b Spin parameter of the halo. 

c Number of SPH particles representing the 
gas. 

d Number of clumps at z — 18.0 for the runs 
marked with "-C" and for all other runs at z — 
37.5. 
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Table 2. List of reactions in our chemical model. 



No. Reaction Rate coefficient (cm 3 s ) Ref. 

1 H + c-^H-+7 h = dcx[- 17.845 + 0.762 log T + 0.1523(logT) 2 1 

- 0.03274(logT) 3 ] T < 6000 K 
= dex[-16.420 + 0.1998(logT) 2 

- 5.447 x 10- 3 (logT) 4 

+ 4.0415 x 10- 5 (logT) 6 ] T> 6000K 

2 H-+H^H 2 + e- fc 2 = 1.5 x T< 300K 2 

= 4.0 x 10- 9 T-°- 17 T > 300 K 

3 H + H+^H++ 7 k 3 = dcx[- 19.38 - 1.523 log T 3 

+ 1.118(logT) 2 - 0.1269(logT) 3 ] 

4 H + H+^H 2 +H+ k 4 = 6.4 x 10~ 10 4 

5 H-+H+^H + H k 5 = 5.7 x lO^T" 05 + 6.3 x 10~ 8 5 

- 9.2 x io- n T - 5 + 4.4 x 10~ 13 T 

6 H++e-^H + H fc 6 = 1.0xl0- 8 T<617K 6 

= 1.32 x io- 6 T-°- 76 T > 617 K 

7 H 2 +H+^H++H k 7 = [-3.3232183 x 10~ 7 7 

+ 3.3735382 x lO^lnT 

- 1.4491368 x 10- 7 (lnT) 2 
+ 3.4172805 x 10- 8 (lnT) 3 

- 4.7813720 x 10- 9 (lnT) 4 
+ 3.9731542 x 10- 10 (lnT) 5 

- 1.8171411 x 10- n (lnT) 6 
+ 3.5311932 x 10- 13 (lnT) 7 ] 
x exp ( ~ 212 T 37 - 15 ) 

8 H 2 +e-^H + H + c- k 8 = 3.73 x l(r 9 T 01121 exp (^o) 8 

9 H 2 +H^H + H + H fcg.i = 6.67 x lO^T 1 / 2 exp [-(1 + ^P)] 9 

fc 9 , h = 3.52 x 10- 9 exp (-«) 10 
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Table 2 — Continued 



No. Reaction 



Rate coefficient (cm 3 s 1 ) 



Rcf. 



3.0 - 0.416 log (i) - 0.327 {log } 2 



n cr ,H = dex 

in h — 5.996xlO~ 30 T'' 1 * / 7.4 \ 
1U il 2 + tl 2 — > Jri 2 + ii + ti fc 10,l — (1,0+6.761 xlO-^T) 8 - 6881 eX P V T J 

fc 10 , h = 1.3 x 10- 9 cxp(-^20) 

n cr , H2 = dex [4.845 - 1.3 log + 1.62 {log (j^)}' 



11 H + c-^H++e-+ 1 



12 D + e-^D++e-+. 



13 H++e-^H + 7 



14 D++e-^D + 7 

15 H~ + c~ —> H + c~ + e" 



fc n = exp[-3. 271396786 x 10 1 
+ 1.35365560 x lOMnTo 

- 5.73932875 x 10°(lnT c ) 2 
+ 1.56315498 x 10°(lnT c ) 3 

- 2.87705600 x lO^lni;) 4 
+ 3.48255977 x 10- 2 (lnT o ) 5 

- 2.63197617 x 10- 3 (lnT c ) 6 
+ 1.11954395 x 10- 4 (lnT o ) 7 

- 2.03914985 x 10- 6 (lnT c ) 8 ] 

k 12 = kn 

k 13 . A = 1.269 x lO-^aifll) 1503 

X [1.0+ (60«25) - 470 ]-1.923 

fc 13 , B = 2.753 x lO" 1 ^ 3 ^) 1 - 500 

X [1.0+ (il5p8) a40 Y 2 - 242 



C14 



415 



«13 



exp[- 1.801849334 x 10 1 
+ 2.36085220 x 10°lnT c 



- 2.82744300 x 10" 1 
+ 1.62331664 x 10~ 2 

- 3.36501203 x 10~ 2 
+ 1.17832978 x 10~ 2 

- 1.65619470 x 10" 3 
+ 1.06827520 x 10~ 4 

- 2.63128581 x 10~ 6 



;inT c ) 2 
;inT c ) 3 
lnT c ) 4 

>T e ) 8 

4nT c ) 6 

>T e ) 7 

;inT c ) 8 ] 



10 

11 

12 
12 
13 



Case A 
Case B 



14 
14 



13 
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Table 2 — Continued 



No. Reaction Rate coefficient (cm 3 s 1 ) Ref. 



16 H-+H^H + H + e- 



17 H" + H+ -> H+ + e- 

18 H + D+^D + H+ 

19 D + H+ -> H + D+ 

20 H 2 + D+ -> HD + H+ 

21 HD + H+ -> H 2 + D+ 

22 H 2 + D -» HD + H 

23 HD + H -» D + H 2 

24 He + e" -> Hc+ + e" + - 



fcie = 2.5634 x 10- 9 T e 178186 
= cxp[-2. 0372609 x 10 1 
+ 1.13944933 x 10°lnT o 

- 1.4210135 x lO-^lnTe) 2 
+ 8.4644554 x 10" 3 (lnT o ) 3 

- 1.4327641 x 10- 3 (lnT c ) 4 
+ 2.0122503 x 10- 4 (lnT c ) 5 
+ 8.6639632 x 10- 5 (lnT c ) 6 

- 2.5850097 x 10- 5 (lnT c ) 7 
+ 2.4555012 x 10- 6 (lnT o ) 8 

- 8.0683825 x 10- 8 (lnT o ) 9 ] 

fci 7 = 6.9 x io- 9 T-° 35 
= 9.6 x l0- 7 T-°- 90 

jfeis = 2.06 x 10- 10 T°- 396 cxp (-f ) 
+ 2.03 x i - 9 T-°- 332 

2.0 x 10- 10 T°- 402 exp(-^i) 

- 3.31 x io-^T 1 - 48 
3.44 x io- 10 T - 35 

[0.417 + 0.846 logT-0.137(logT) 2 ] x 10~ 9 

1.1 x 10- 9 cxp(-^) 



T P < 0.1 eV 



13 



kig = 

k 2 o = 
&21 

k 2 2 



k 



23 



4680 



1.69 x 10- 10 exp(- 
1.69xl0- 10 exp(-^2 + i^) 

5.25 x 10- n exp(-^f2) 

5.25 x lO-Hexpf-^ + i 7 ^) 



fc 24 = cxp[-4.409864886 x 10 1 
+ 2.391596563 x lO^nTc 
- 1.07532302 x lO^lni;) 2 



T c > 0.1 eV 

T < 8000 K 
T > 8000 K 



T < 2 x 10 5 K 
T > 2 x 10 5 K 



T < 200 K 
T > 200 K 

T < 200 K 
T > 200 K 



15 
16 
16 

17 
17 
18 

19 

13 
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Table 2 — Continued 



No. Reaction 



Rate coefficient (cm 3 s x ) 



+ 3.05803875 x 10°(lnT o ) 3 

- 5.6851189 x lO-^lnTo) 4 
+ 6.79539123 x 10- 2 (lnT o ) 5 

- 5.0090561 x 10- 3 (lnT o ) 6 
+ 2.06723616 x 10- 4 (lnT c ) 7 

- 3.64916141 x 10- 6 (h N 

25 He+ + c- -» He + 7 fc 25 ,rr,A = lO" 11 ?" - 5 



Rcf. 



T e f] 



10 -n T -o.5 [ 12 _ 72 _ 1.615 logT 
- 0.3162(logT) 2 + 0.0493(logT) 3 ] 



fc25,rr,B = lO" 1 ^" - 5 [11.19 - 1.676 log T 

- 0.2852(logT) 2 + 0.04433(logT) 3 ] 



fc 25 di = 1.9 x 10- 3 T- 15 exp 
x [1.0 + 0.3 exp 



31 Si+ + c" 



Si + 7 



= 1.23 xl0-^(^) 2 49 exp (^p) 
= 9.62 x 10- 8 (^^cxp ( - 115 T 786 - 2 

fc 31 =7.5xl0- 12 (^)- ' 55 
= 4.86xl0--(^)-°- 32 
= 9.08xl0--(^) ' 818 



) 



Case A 



Case B 



20 



20 



21 



26 


He+ + H - 


-> He + H+ 


he 


= 1.25xl0--(^)°' 25 




22 


27 


He + H+ - 


-» Hc + + H 


&27 


= 1.26 x 10-9T- - 75 exp ( 127 T 500 ) 
= 4.0 x i(r 37 T 4 - 74 


T < 10000 K 
T > 10000 K 


23 


28 


He+ + D - 


-» He + D+ 


hs 


= he 






29 


He + D+ - 


-> He + + D 




= hi 






30 


C+ + o- - 


♦ C + 7 


ho 


= 4.67xl0--(^)-°- 6 


T < 7950 K 


24 



7950 K < T < 21140 K 
T > 21140 K 

T < 2000 K 

2000 K < T < 10 4 K 

T > 10 4 K 



25 



32 0+ + < 



+ 7 



V32 



1.30 x io- 10 T-°- 64 

1.41 x io- 10 T-°- 66 + 7.4 x 10- 4 T- 15 
x exp (-^m) [1.0 + 0.062 x exp (_1«QQ0 



T < 400 K 
_)] 7 101) K 



26 
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Table 2 — Continued 



No. Reaction 



Rate coefficient (cm 3 s 1 ) 



Rcf. 



33 C +e" -» C+ +e" + e" 

34 Si + e" -> Si+ + e" + e" 

35 O + c- -> 0+ +e- + e" 

36 0+ + H -» O + H+ 

37 O + H+ -» 0+ + H 



fc 33 = 6.85 x 10" 8 (0.193 + u)- 1 ^ - 2 ^-" 



u = 11.26/7; 



fc 34 = 1.88 x 10" 7 (1.0 + u°- 5 )(0.376 + u)-^ - 2 ^-" u = 8.2/T e 



&35 
k 3 6 



3.59 x 10- 8 (0.073 + -u)- 1 ?i - 34 e- u 

4.99 x io- n T - 405 + 7.54 x lO" 10 ! 1 " - 458 



k 37 = [1.08 x lo-nr - 517 

+ 4.00 x iQ-iOT - 00669 ] cxp (-^?) 



38 + Hc+^0+ + He k 38 = 4.991 x lO" 15 (n^oo)"' 3794 cxp (-^ 



21000 J 



39 C + H+ -» C+ + H 

40 C+ + H -> C + H+ 

41 C + Hc+ -» C+ + He 

42 Si + H+ -> Si+ + H 

43 Si + He+ -» Si+ + He 

44 C+ + Si -» C + Si+ 



k 3 9 = 



+ 2.780 xl0" 15 (i)- ■ cxp( M ^ M ) 
3.9 x io-^t - 213 



k 40 = 6.08 x 10- 14 (t^) 1 - 96 exp (-™) 



fc 4 i = 8.58 x i(T 17 T - 757 
= 3.25 x i(T 17 T -9 68 
= 2.77 x io- 19 T 1597 

k i2 = 5.88 x 10-13T - 848 
= 1.45 x 10~ 13 T 

k i3 = 3.3 x lO- 9 

k i4 = 2.1 x 10~ 9 



45 Si+ + H+ -» Si++ + H k 45 = 4.10 x HT 10 (j^) 



0.24 



1.0 + 3.176^(53^)] exp (-^) 



46 Si++ + H -» Si+ + H+ fc 46 - 1.23 x lO" 9 

x 



47 Si++ + c- -» Si+ + 7 



1.0 + 3.17 exp (^xw) 

, ,„ / T \ -0.6346 

/c 47 , rr = 1-75 x 10" 12 (-^ ) 



u = 13.6/Te 



T < 200 K 

200 < T < 2000 K 

T > 2000 K 

T < 10 4 K 
T > 10 4 K 



27 
27 
27 
28 
29 

30 

29 
29 
31 

32 

33 
33 
32 

32 

34 
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Table 2 — Continued 

No. Reaction Rate coefficient (cm 3 s -1 ) Ref. 

fc 47 ,di = 2.2552 x lO" 11 ^" 1 - 5 exp (-^p) 35 
+ 5.6058 x 10- 9 T- 15 exp (-^) 

75 H 2 + Hc^H + H + He fe 75| i = dex [ — 27.029 + 3.801 log (T) — 29487/T] 36 

/c 75ih = dcx [-2.729 - 1.75 log (T) - 23474/T] 

n crHc = dex [5.0792(1.0 - 1.23 x 10~ 5 (T - 2000)] 36 

76 OH + H^O + H + H fc 76 = 6.0 x lO" 9 exp (-^o) 33 

77 HOC+ + H 2 -» HCO+ + H 2 fc 77 = 3.8 x lO" 10 37 

78 HOC+ + CO -» HCO+ + CO fc 78 = 4.0 x lO" 10 38 

79 C + H 2 ^CH + H fc 79 =6.64 x 10- 10 exp (-i^o) 39 

80 CH + H^C + H 2 fc 80 = 1.31 x 10- 10 exp (-f) 40 

81 CH + H 2 ^CH 2 +H fc 8 i =5.46 x 10- 10 exp 41 

82 CH + C^C 2 +H fc 82 = 6.59 x 10~ n 42 

83 CH + O^CO + H fc 83 = 6.6 x lO" 11 T< 2000K 43 

= 1.02 x 10~ 10 exp (-^) T> 2000K 44 

84 CH 2 + H -> CH + H 2 fc 84 = 6.64 x lO" 11 45 

85 CH 2 + O^CO + H + H fc 85 = 1.33 x 10~ 10 46 

86 CH 2 + O^CO + H 2 fc 86 = 8.0 x lO" 11 47 

87 C 2 + O^CO + C fc 87 = 5.0 x 10- 11 (^q)°' 5 T< 300K 48 

= 5.0 x 10- 11 (^o) ' 757 T > 300 K 49 

88 + H 2 ^OH + H fc 88 = 3.14 x 1Q- 13 (-^) 2 ' 7 exp (-^) 50 
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Table 2 — Continued 

No. Reaction Rate coefficient (cm 3 s _1 ) Ref. 

89 OH + H^O + H 2 /c 89 = 6.99x 10- 14 (3^) 2,8 cxp(-if2) 51 

90 OH + H 2 ^H 2 + H /c 90 = 2.05 x 10- 12 (3^) 152 exp(-if^) 52 

91 OH + C^CO + H jfcgi = 1.0 x 10- 10 42 

92 OH + 0^0 2 + H k 92 = 3.50 x 10~ n T < 261 K 53 

= 1.77 x 10- n exp(^) T>261K 41 

93 OH + OH^H 2 + H k 93 = 1.65 x 10~ 12 (-^) 1,14 exp (-f ) 42 

94 H 2 + H^H 2 + OH k 94 = 1.59 x lO" 11 (^) 1-2 exp (-^f 2 ) 54 

95 2 + H^OH + /c 95 = 2.61 x 10- 10 cxp(-^) 41 

96 2 + H 2 ^OH + OH /j 96 = 3.16 x 10- 10 exp (-^0) 55 

97 2 + C^CO + /c 97 = 4.7x lO- 1 ^^) -0 ' 34 T< 295K 42 

= 2.48 x 10~ 12 (ggjj) 1,54 exp T > 295 K 41 

98 CO + H^C + OH /fc 98 = 1.1 x 10- 10 (^)°- 5 exp (-1^00) 33 

99 H++H 2 ^H++H fc 99 = 2.24 x 10- 9 (3^)° ° 42 exp (-3^) 56 

100 H++H^H++H 2 fc 100 ^ 7.7 x 10- 9 cxp(-i^o) 57 

101 C + H+^CH++H fcioi = 2.4 x 10~ 9 33 

102 C + H+ -» CH+ + H 2 fci 02 = 2.0 x 10~ 9 33 

103 C++H 2 ^CH++H fci 03 = 1.0 x lO" 10 exp (-^jP) 58 

104 CH++H^C++H 2 k 104 = 7.5 x lO" 10 59 

105 CH+ + H 2 -> CH+ + H fcios = 1.2 x 10~ 9 59 

106 CH+ + O^CO+ + H fcioe = 3.5 x lO" 10 60 
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Table 2 — Continued 

No. Reaction Rate coefficient (cm 3 s _1 ) Ref. 

107 CH 2 + H+ -> CH+ + H 2 fc 107 = 1.4 x lO" 9 33 

108 CH++H^CH++H 2 fcios = 1.0 x 10" 9 exp(-^) 33 

109 CRt + H 2 -» CHj + H fciog = 1.6 x lO" 9 61 

110 CHj + O -» HCO+ + H fcno = 7.5 x 10~ 10 33 

111 CH++H^CH++H 2 fcm = 7.0 x 10- 10 exp(-i^) 33 

112 CH+ + O -» HCO+ + H 2 fc n2 = 4.0 x 10- 10 62 

113 C 2 + 0+^CO+ + C fcn3 = 4.8 x 10~ 10 33 

114 0++H 2 ^OH++H fcn4 = 1.7 x 10- 9 63 

115 + H+^OH++H fcii 5 = 1.5 x 10" 9 33 

116 + H+^OH++H 2 fcne = 8.4 x 10- 10 64 

117 OH + H+ -» H 2 0+ + H 2 fciiy = 1.3 x 10" 9 33 

118 OH + C+^CO++H fc n8 = 7.7 x 10- 10 33 

119 OH+ +H 2 H 2 0+ +H fcng = 1.01 x 10~ 9 65 

120 H 2 0+ + H 2 -» H 3 0+ + H fci 20 = 6.4 x 10- 10 66 

121 H 2 + H+ -> H 3 0+ +H 2 fci2i = 5.9 x 10~ 9 67 

122 H 2 + C+ ^ HCO+ +H fc 122 = 9.0 x lO" 10 68 

123 H 2 + C+ -» HOC+ +H fci 23 = 1.8 x 10~ 9 68 

124 H 3 + + C^HCO++H 2 fc m - 1.0 x lO" 11 33 
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Table 2 — Continued 



No. Reaction 

125 H 3 0+ + Si -» H 2 + SiH+ 

126 2 + C+ -» CO+ + O 

127 2 + C+ -» CO + 0+ 

128 2 + CH+ -» HCO+ + OH 

129 0+ + C -» C0+ + O 

130 CO + H+ -» HOC+ + H 2 

131 CO + H+ -» HCO+ + H 2 

132 HCO+ + C -» CO + CH+ 

133 HCO+ + Si -» CO + SiH+ 

134 HCO+ +H 2 -» CO + H3O+ 

135 Si + H+ -» SiH+ + H 2 

136 SiH+ + H -» Si+ + H 2 

137 SiH+ +H 2 -» Si + H 3 0+ 

138 H 2 + Hc+ -» He + H+ 

139 H 2 + He+ -» He + H + H+ 

140 CH + H+ -> CH+ + H 

141 CH 2 + H+ -> CH+ + H 

142 CH 2 + He+ -» C+ + He + H 2 

143 C 2 + He+ -» C+ + C + He 



Rate coefficient (cm 3 s 1 ) Ref. 

fci25 = 1.8 x 10~ 9 33 

ki2& = 3.8 x 10- 10 61 

fci 27 = 6.2 x 10- 10 61 

ki2s = 9.1 x 10- 10 61 

fci 29 = 5.2 x 10~ n 33 

fci 30 = 2.7 x 10- 11 69 

fci3i = 1-7 x 10- 9 69 

fci 32 = 1.1 x 10- 9 33 

fciss = 1-6 x 10~ 9 33 

fci34 = 2.5 x 10~ 9 70 

ki35 = 3.7 x 10 - 9 33 

fci 36 = 1-9 x 10~ 9 71 

fc m = 8.0 x 10- 10 71 

fci 38 = 7.2 x 10~ 15 72 

fc 139 =3.7x 10- 14 cxp(f ) 72 

k 140 = 1.9 x 10- 9 33 

fcui = 1.4 x 10- 9 33 

k li2 = 7.5 x 10~ 10 33 

feus = 1.6 x 10~ 9 33 
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Table 2 — Continued 

No. Reaction Rate coefficient (cm 3 s _1 ) Ref. 

144 OH + H+^OH++H fc 144 = 2.1 x 1CT 9 33 

145 OH + Hc+ -» 0+ +He + H k li5 = 1.1 x 10" 9 33 

146 H 2 + H+ -» H 2 0+ +H fc 146 = 6.9 x 1(T 9 73 

147 H 2 + Hc+ -» OH + Hc + H+ fc 147 - 2.04 x KT 10 74 

148 H 2 + Hc+ -» OH+ +He + H fc 148 = 2.86 x KT 10 74 

149 H 2 + He+ -» H 2 0+ + He fci 49 = 6.05 x lO" 11 74 

150 2 + H+ -> + H fciso = 2.0 x lO" 9 73 

151 2 + Hc+ -> + He fci 5 i - 3.3 x lO" 11 75 

152 2 + Hc+^0+ + + Hc fci 52 = 1.1 x 10~ 9 75 

153 0^ + 0^02 + 0+ fci 53 = 5.2 x 10~ n 33 

154 CO + Hc+ -» C+ +0 + Hc k 154 = 1.4 x 10~ 9 (^)~°' 5 76 

155 CO + Hc+ -> C + 0+ +Hc fc 155 = 1.4 x 10" 16 (^)~°' 5 76 

156 CO++H^CO + H+ fci 56 = 7.5 x 10~ 10 77 

157 C- + H+ -> C + H k 157 = 2.3 x lO" 7 (^) "°' 5 33 

158 CT + H+ -» O + H fci 58 = 2.3 x 10~ 7 (^) "°' 5 33 

159 He+ + H-^Hc + H fc 159 = 2.32 x lO" 7 (^y)~°' 52 exp 78 

160 H+ + e" -» H 2 + H fci 60 = 2.34 x 10~ 8 (^>o)~ ' 52 79 

161 H++e-^H + H + H k W i = 4.36 x lO" 8 (^>o)~ ' 52 79 
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Table 2 — Continued 



No. Reaction Rate coefficient (cm 3 s 1 ) Ref. 

162 CH++c-^C + H A:i62 = 7.0 x 10- 8 (g^) -0 ' 5 80 

163 CH++c-^CH + H jfci 63 = 1.6 x lO" 7 (g^) -0 ' 6 81 

164 CH++c-^C + H + H /ci64 = 4.03x 10- 7 (3^) _0 ' 6 81 

165 CH++c-^C + H 2 fci 65 = 7.68 x 10~ 8 (^) ~°' 6 81 

166 CH++c-^CH 2 +H fci 66 = 7.75 x 10~ 8 (^) "°' 5 82 

167 CH++c-^CH + H 2 fc 167 = 1.95 x lO" 7 ^)" 05 82 

168 CH+ + e" -» CH + H + H k 16R = 2.0 x 10~ 7 (^o)~°' 4 33 

169 OH++c-^0 + H /ci 69 = 6.3 x 10- 9 (g^) -048 83 

170 H 2 0+ + c- -> O + H + H /ci 70 = 3.05 x 10- 7 (^q)" ' 5 84 

171 H 2 + +e-^0 + H 2 k 171 = 3.9 x 10~ 8 (^o)~°' 5 84 

172 H 2 0+ + e- ^ OH + H k 172 = 8.6 x 10~ 8 (^)~°' 5 84 

173 H 3 0+ + c" -> H + H 2 /ci7 3 - 1.08 x 10- 7 (g^) -0 ' 5 85 

174 H 3 0+ + c-^OH + H 2 fci 74 = 6.02 x 10 ~ 8 (^) "°' 5 85 

175 H 3 + +e" -» OH + H + H fc 175 = 2.58 x 10 ~ 7 (^)~°' 5 85 

176 H 3 + +e~ -> + H + H 2 fci 76 = 5.6 x 10~ 9 (^o)~°' 5 85 

177 O^+c-^O + O fcir? = 1-95 x 10~ 7 (^) 86 

178 CO++c-^C + fci78 = 2.75x 10- 7 (3^) -0 - 55 87 

179 HCO+ + e~ -» CO + H fci 79 = 2.76 x 10~ 7 ( ^) ~°' 64 88 

180 HCO++e-^OH + C fci 80 = 2.4x lQ-^^)^ 64 88 
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Table 2 — Continued 

No. Reaction Rate coefficient (cm 3 s _1 ) Ref. 

181 HOC+ + c-^CO + H fcisi = 1.1 x lO" 7 ^) -10 33 

182 SiH+ + c" -» Si + H k 1S2 = 2.0 x 1(T 7 (^) _ °' 5 33 

183 H-+C^CH + e- fci 83 = 1.0 x 1(T 9 33 

184 H-+O^OH + e- fcis4 = 1.0 x 10~ 9 33 

185 H-+OH^H 2 + e" fci 85 = 1.0 x 10" 10 33 

186 CT+H^CH + e" fc 186 = 5.0 x lO" 10 33 

187 C-+H 2 ^CH 2 +c- fc 187 = 1.0 x 10~ 13 33 

188 CT+O^CO + e" fci 88 = 5.0 x lO" 10 33 

189 CT+H^OH + e- fc 189 = 5.0 x 10~ 10 33 

190 O" +H 2 ->H 2 + e- k 190 = 7.0 x lO" 10 33 

191 CT+C^CO + e" feigi = 5.0 x 10- 10 33 

192 H 2 +H+^H++7 k 192 = 1.0 x 10~ 16 89 

193 C + c-^C-+ 7 fci93 = 2.25 x 10~ 15 90 

194 C + H^CH + 7 fci 94 = 1.0 x 10~ 17 91 

195 C + H 2 ^CH 2 + 7 fci 95 = 1.0 x 10" 17 91 

196 C + C^C 2 + 7 fc 196 = 4.36 x 10- 18 (^)°' 35 cxp (-if 3 ) 92 

197 C + O^CO + 7 fci 97 = 2.1 x 10~ 19 T< 300K 93 

= 3.09x 10- 17 (3^)°- 33 exp(-lf^) T> 300K 94 

198 C++H^CH++7 fc 198 = 4.46 x 10- 16 T-°- 5 exp(-|||) 95 
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Table 2 — Continued 

No. Reaction Rate coefficient (cm 3 s _1 ) Ref. 

199 C++H 2 ^CH++ 7 fci 99 =4.0x 10- 16 (3^)" 02 96 

200 C+ + O^CO++7 fc 200 = 2.5 x 10" 18 T< 300K 93 

= 3.14x 10- 18 (^ J )" 015 cxp(f ) T> 300K 

201 + c-^0-+7 fc 20 i = 1.5 x 10~ 15 33 

202 O + H^OH + 7 fc 202 = 9.9 x 10- 19 (3^) _ ° 38 33 

203 + ^0 2 +7 fc 2 o 3 =4.9x 10- 20 (3^) 158 91 

204 OH + H^H 2 + 7 fc 204 = 5.26 x lO" 18 (^jo)~ 5 ' 22 exp (-f ) 97 

205 Si++H^SiH++ 7 k 205 = 1.17 x 10~ 17 (^o)~°' 14 98 

206 H + H + H^H 2 +H fc 206 = 1.32 x lO" 32 (^o)~°' 38 T< 300K 99 

= 1.32 x 10" 32 (^o)" 10 T> 300K 100 

207 H + H + H 2 ^H 2 + H 2 fc 207 = 2.8 x lO" 31 ! 1 - - 6 101 

208 H + H + He -> H 2 + He fc 208 = 6.9 x 10- 32 T-° A 102 

209 C + C + M^C 2 + M fc 209 = 5.99 x 10 ~ 33 (glL) -1-6 T< 5000K 103 

= 5.99 x HT 33 (^)- - 64 exp (^») T> 5000K 104 

210 C + + M->CO + M k 210 = 6.16 x lO" 29 (^y) _3 '° 8 T< 2000K 105 

= 2.14x 10- 29 (^ ?J )^ 3 08 exp( 2 ^) T> 2000K 76 

211 C+ + + M -> CO+ + M k 2n = 100 x k 2W 76 

212 C + 0+ + M -> CO + + M fc 2 i 2 = 100 x k 2W 76 

213 + H + M->OH + M k 213 = 4.33 x lO" 32 51 

214 OH + H + M ^ H 2 + M k 21i = 2.56 x 10~ 31 (^o)" 2 '° 105 



-39- 



Table 2 — Continued 



No. Reaction 



Rate coefficient (cm 3 s 



3 a -l\ 



Ref. 



215 + + M-> 2 + M fc 2 i 5 = 9.2xl0- 



^300/ 



45 



216 + CH^HCO++e" fc 216 = 2.0 x lO" 11 (^)°' 44 106 
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Table 3. List of processes in our thermal model. 



Process 



References 



Cooling: 

Lyman-o cooling 

He electronic excitation 

Thermal bremsstrahlung 

Compton cooling 

H2 rotational, vibrational lines 

HD rotational, vibrational lines 

Fine structure lines (C, C + , O, Si, Si + ] 

CO rotational, vibrational lines 

H2O rotational, vibrational lines 

OH rotational, vibrational lines lines 

H + recombination 

He + recombination 

H k He collisional ionization 

H2 collisional dissociation 

Heating: 

H2 gas-phase formation 



Cen (1992) 
Cen (1992); Brav et al. (2000) 
Shapiro k Kang (1987 ) 

Cen (1992) 
Le Bourlot et al. (1999) 
Lipovka. Nunez-Lopez, k Avila-Reese (2005) 
Many sources; see Paper I 
Neufeld k Kaufman (1993); Neufeld. Lepp k Melnick (1995. ) 
Neufeld k Kaufman (1993); Neufeld. Lepp k Melnick (1995) 
Pavlovski et al. (2002) 
Ferland et al. (1992) 
Aldrovandi k Peauignot (1973 ); Hummer k Storey (1998) 
Janev et al. (1987) 
See g2] 

Many sources; see Paper I and $2] 



